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1 Introduction 



The first version of tlie Pair Boson Hamiltonian (PBH) model was proposed by Zubarev 
and Tserkovnikov in 1958 p. Tlieir intention was to generalize the Bogoliubov model of the 
Weakly Imperfect Bose Gas j2I by including more terms from the total interaction, without 
losing the possibility of having an exact solution. We refer the reader to [3] and to |3] for a 
more recent discussion of this question. 

The suggestion of Zubarev and Tserkovnikov was to consider a truncated Hamiltonian 
which includes a diagonal term representing forward-scattering and exchange-scattering as 
well as a non-diagonal BCS-type interaction term. The model containing only the forward- 
scattering part of the interaction corresponds to the Mean-Field (or the Imperfect) Bose gas, 
see jl] and j3] for details. Using the same method as they had used earlier for the fermion 
BCS model jH], the authors give in pP a "solution" of the PBH model. Later this Hamil- 
tonian became the subject of very intensive analysis (Zj-jH], leading essentially to the same 
conclusion as in P , namely, that the PBH has the same thermodynamic properties as a cer- 
tain approximating Hamiltonian quadratic in the creation and annihilation operators. Using 
this Hamiltonian which can be diagonalized by the canonical Bogoliubov transformation, 
its thermodynamic properties were investigated and it was shown to have some intriguing 
properties. One of these is possibility of the occurrence of two kinds of condensation, the 
standard one-particle Bose-Einstein condensation as well as a BCS-type pair condensation 
which may appear in two stages, see e.g. fOj, Another one concerns the gap in the 

spectrum of "elementary excitations" [Z]-|n]- In spite of fairly convincing arguments these 
papers did not prove rigorously that the above mentioned solution of the PBH model is exact. 
A mathematical treatment of the PBH model, related to representations of the Canonical 
Commutation Relations (CCR) appeared in |12j . 

In the present paper we give a variational formula for the pressure for the PBH model and 
provide a rigorous derivation of the formula. The latter yields the same expression for the 
pressure as was obtained in the corresponding Euler-Lagrange equations coinciding with 
self-consistency equations studied in P and |Z|-|I2I- In an earlier paper jTSj we conjectured 
that the pressure can be expressed as the supremum of a variational functional depending 
on two measures: a positive measure describing the particle density and a complex measure 
describing the pair density, similar to the Cooper pairs density in the BCS model. This con- 
firmed the conclusion of ^U], about the coexistence of one-particle and pair condensates. 
The study in ^Hj was inspired by the Large Deviation Principle (LDP) developed for the 
analysis of boson systems in This method gives rigorous results for the pressure 

in the case of models with diagonal (commutative) boson interactions. A similar technique 
was developed in |THj-[221 based on the work [22j, extending the LDP to noncommutative 
Mean-Field models (including the BCS one) with only bounded operators involved in Hamil- 
tonians. Since neither of these methods apply to the PBH without extensive modifications, 
here we opted for the Approximating Hamiltonian Method (AHM) [25; which has been al- 
ready successfully applied to many models, including some interacting boson models (see for 
example U, 0, |2S1). 

There is renewed interest in the properties of the PBH interaction in the context of finite 
boson systems confined in a magneto-optic trap, see e.g. [2n]-|2H]- We do not discuss this 
aspect in the framework of our approach leaving it for future publications. 

Now we turn to the exact formulation of the PBH model in its simplest form, that is, with 



constant pair and mean- field boson couplings [T!^ . 



Let A C M*^ be a cube of volume = |A| centered at the origin. Then the kinetic energy oper- 
ator for a particle of mass m confined to the cubic box A, that is the operator — A/2m with pe- 
riodic boundary conditions, has eigenvalues e{k) = \\k\\'^/2m, k G A* := {2its/V^^'^\s E Tj"}. 
Consider a system of identical bosons of mass m enclosed in A. For k E A* let and a^, be 
the usual boson creation and annihilation operators satisfying the CCR [a^ , a^/] = 5^ ^/ and 
let Nk := a^a^ be the fc-mode particle number operator. The kinetic-energy operator Ta for 
the Perfect Bose-gas, can be expressed in the form := YlkeA* 
To introduce a pairing term in the Hamiltonian we shall need the operators 

Ak = A_k ■■= a^a.k , k e A* . (1.1) 

Let 

iVA := Yl •= Yl ^(^)^'^ ' (1-2) 

fceA* fceA* 

where the function A : M'' i— C satisfies the following conditions: 

|A(A;)| < |A(0)| = 1 , ~X{k) = ~X{-k) for all keW , 
there exists C < oo and 6 > such that 

l-^^^^l - I + ||^||max(i.,i./2+l)+6 (^-2) 



for all k e W. Note that (fOjl implies that A G L^(W) and that there exists M < oo such 
that 

m^- J2 \kk)\<MV, (1.4) 



n^:= Y,<k)\\{k)\^ <MV, (1.5) 



and 



fceA* 



snp e{k)\X{k)\^ < M (1.6) 
fceA* 



for all AcW. 

Then for constant couplings u, v the PBH is defined by 



H.:=T,-^Q:Q. + ^NI (1.7) 



Remark 1.1 Let ip := argA(O) and X{k) := X{k)e Then A(0) = 1 and we can write Hf, 
in the form 

H. = T,-—QIQ, + —Nl (1.8) 

with 

Q,:=YKk)A,, (1.9) 

fceA* 

where \X{k)\ < A(0) = 1 for all k e W. 



Remark 1.2 We shall assume that v > and a := v — u > 0. The latter condition ensures 
the superstability of the model, see Theorem \2. 1\ Note that in the case m < {BCS repulsion) , 
the second condition a > is trivially satisfied. In we have proved that the case u < 
gives the same thermodynamics as the Mean-Field (MF) Bose-gas: 

Hr--=T. + ^Nl. (1.10) 

Thus in deriving the variational formula we emphasize the case m > 0. We recall that this 
condition is necessary for nontrivial condensation of boson pairs, see e.g. J^-USf. We shall 
discuss the relation between these conditions and the thermodynamic properties of the model 
M.S^} in Section\^ 



For the convenience of the reader we now state (without proof) the principal theorems and 
describe the the logical sequence used in proving the main result of this paper. We shall 
need the grand-canonical pressures for several approximating Hamiltonians. Recall that for 
an inverse temperature (3 and a chemical potential /i the the grand-canonical pressure for a 
system with Hamiltonian Ti 

A IS 

l-lnTiexp{-f]{n^-fiN^)} . (1.11) 

For simplicity in the sequel we shall omit the thermodynamic variables (3 and fi and we shall 
write, for example, Pa for the grand-canonical pressure corresponding to the Hamiltonians 

:= ^ InTi exp {-P{H^ - fiN^)}. (1.12) 
We shall denote the thermodynamic limit A | M*^ by the symbol ' lim '. 

A 

Consider the approximating Hamiltonian 

i/f (g, p) := Ta + vpN, - ^u{Qlq + Q,q*) - ^ V + ju\q\' , (1.13) 

where g G C and p G M+ are variational parameters. The Hamiltonian H^^ (g, p) can be 
diagonalized and the corresponding pressure p\ (g, p) can be calculated explicitly to give in 
the thermodynamic limit 

p^^\q,p): = limpi^^(g,p) 

A 

L 1 in[l _ exp{-(3E{k,q,p))] - ^ {E{k,q,p) - f{k,p)) 



where 
with 



E{k,q,p):={f\k,p)-\h{k,qW/' , (1.15) 



f{k,p):=e{k)—p + vp and h{k,q) := uq \*{k) . (1-16) 
Using (jl.l3p the Hamiltonian ()1.8|) can be written identically as 

H, = Hf\q,p) + Hl{q,p) (1.17) 



where 



Hl{q. p) ■■= -^<Ql - y^HQ. - Vq) + - pf . (1.18) 



The main result of this paper states that if the variational parameters q and p are chosen 
in an "optimal" way, then the contribution to the pressure arising from the residual term 
H^{q,p) vanishes in the thermodynamic hmit. 

Let us define the following function for q > and p > 

a(g,p):= inf {f{k,p)-\h{k,q)\)=vp-p-\u\q, (1.19) 



see 



dnni). 



Theorem 1.1 The limiting pressure for the PBH model il.^) with u > {BCS attraction) 
has the form 

p := limpA = sup inf p*-^-* (g, p) = sup inf p^'^\q,p) , (1-20) 
while with u < {BCS repulsion) it has the form 

p := limpA = inf inf p''^^(g, p) = inf inf p^'^\q, p) ■ (1-21) 

A (jGCp>0 '3>0p;o-(<j,p)>0 

Note that to obtain the approximating Hamiltonian (jl.l3p . the term —uQ*jjQ^/2V in p.8p 
is replaced by -u{Qlq + Q^q*)/2 + Vu\q\'^/2 and vNl/2V by vpNj, - Vvp^/2. 

We shall prove Theorem 11.11 in two steps. Here we describe these steps for m > and before 
the end of the section we indicate the modifications necessary for the case u < 0. 

The first step which we call the first approximation is to linearize the term —uQ\Q J2'V in 
i^A- For technical reasons we need to add to our Hamiltonians some source terms. Therefore, 
we define for z/, 77 G C 



H^{iy, ri) ■= - {uQl + z/*g J - ^/V {r^al + r^*a^) , (1.22) 
and the first approximating Hamiltonian 

H^\q, V, ri) := + ^iV^ - \u{Q\q + Q,q*) + ]A/u\qf - (1.23) 



{vQ\ + z/*g J - (r/aS + v^a^) . 

From (fr22|l and (fr28|l we have 

H,{v,ri) = H'f^{q,v,ri)^El{q) 

where 

ii'M) = -^^(Ql - VqUQa - Vq) < 0. (1.24) 

First we show (see Section E)) that with the right choice of the parameter q = q, the residual 
perturbation H2{q) does not contribute to p^{i','r]), the pressure for the PBH p.22|) in the 
thermodynamic limit, i.e., the pressure corresponding to the Hamiltonian Ht^{v,ri) coincides 
with the limit of pT* {q-, v)j the pressure for H^\q, u, rj): 



Theorem 1.2 For any v and rj with \u\ < 1 and \r]\ < 1, 



In particular 



\impJu,ri) = limsuppA Vg, i^, 77). (1-25) 

A A g 



limp J rj) = Mm sup p^l\q,ri). (1.26) 

A A q 



where Pa{v) •= Pa{0,v) and p^^\q,r]) := P^a\<1,0,v) ^^^^ the pressures corresponding to the 
Hamiltonians HA{ri) := HA{0,ri) and Hl^\q,ri) := ifi^''(g, 0, 77) respectively. 

Next, in Section |3] we study a second approximating Hamiltonian obtained from ()1.23|) by 
replacing the term vN^/2V by a linear term vpNjy — Vvp'^/2: 

H^^\q, p, r/) := + vpN^ - ^uiQlq + Q^q*) - jvp' + ju\q\^ - (r/a* + . (1.27) 

We denote the pressure corresponding to the Hamiltonian ()1.27p by p^^^ (g, p, rj) . Note that 
by frT!!|l and fTTTjl one has H'f'{q,p,{)) = Hj^\q,p). We shall show in Lemma lO that 
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Pa {q,p,v) =Pa {q,p) + |?7l 

where 6 := argg and ip := argr/. 

Our next theorem establishes a similar variational relation between the pressure Pa{v) ^"^^ 

PaHq^p^v)- 



Theorem 1.3 



limpA('7) = limsup inf Pa ^(g, p, 77) = lim sup inf pi^-* (g, p, 77) , (1.28) 

A A p>0 A g>o p>0 



where for q > we put 

p^I\q,p,v) :=pf (ge^(-+^^),p,r/) = p^^\q, p) + ^III . (1.29) 

/(0,p) - uq 

Note that the difference between the statement in Theorem 11.11 and that in Theorem II. HI 
(apart from the 77 dependence) is that the thermodynamic limit is taken after taking the 
infimum over p and the supremum over g. In the next theorem we show that the order of 
the thermodynamic limit and taking the infimum and supremum can be reversed: 

Theorem 1.4 For rj ^ 0, 

p{ri) := limp A{ri) = snp inf P^'^\q,p,v): (1.30) 

A q>0 P-o-{q,p)>0 



where we put 



^ f{0,P)-uq 



cf. expression il.29\) . 



In Lemma (4. 51 we prove that p = lim^_»op(?7) so that Theorem II .41 gives 

p = hmsup inf p^'^\q, p,r]). (1.32) 

V-*0 q>0 P-o-iq,p)>0 

Finally in Lemma 14.61 we prove that the order of the limit 77 ^ and taking the infimum 
and supremum can be reversed to yield the main result Theorem II. II for the BCS attraction. 

The important difference for the repulsive case, m < 0, is that instead of p.24j) we now have 



Therefore the first approximation (Section IH)) should be constructed in the same way as the 
second approximation (Section 0]). The proof of the second part of Theorem II. H ()1.2H1 . for 
M < is given in Section |S1 (f). 

It is important to note that the variational formula conjectured in jOj has the same Euler- 
Lagrange equations as those given by Theorem 11.11 Thus the detailed study of these equa- 
tions carried out in ^Hj applies to our result. In particular, this concerns the sequence 
of phase transitions in the PBH model ()1.8|) and the conditions for the coexistence of the 
generalized Bose condensation and the condensation of boson pairs, see also Sectional 

The paper is organized as follows. We start by proving in Section |21 that the PBH model 
fll.8|l is superstable. In Sections El and E] we shall assume that u > 0. Section IHl is devoted to 
establishing the first approximation giving the proof of Theorem 11.21 In Section |3] we turn 
to the second approximation giving the proof of Theorem 11.31 and the other results needed 
to obtain Theorem 11.11 for m > 0. Finally in Section El we discuss the variational problem as 
well as related open questions for all values of u and we finish the proof of Theorem 11.11 for 
n < 0. Some commutator relations are given in Appendix A and in Appendix B we give a 
bound needed in our proofs. 

2 Superstability 

In this section we establish the superstability of the PBH model ()1.8j) . When m < super- 
stability is obvious. To prove it for m > and a = v — u > 0, we shall need the following 
lemma which is used in several other places in the paper. 

Lemma 2.1 The following inequality is satisfied 



Lu(Ql-Vq*){Q,-Vq)>0. 



(1.33) 



QIQ^ < Nl + MVN^. 



(2.1) 



Proof: The inequalities 



{X*{k)ak'al ± X*{k')a*_f^,a_f?)* {X*{k)ak'al ± A*(A:')«-fc'«-fc) > 



and definition (jl.lj) imply that for k {k' , —k'}, 

-{Nk + \m\)m> - {N.k' + \Kk')\)N.k 

< -\X{k)\\Nk + l)Nk' - \X{k')\\N^k' + l)N_k 
< X*{k)X{k')AlAk' + X*{k')X{k)Al,Ak 

< \X{k)\\Nk + l)Nk' + \X{k')\^{N.k' + l)iV_fc 
< {Nk + \Xik)\)Nk, + iN_k, + \X{k')\)N_k . 



(2.2) 



By P-lj) we also have 



AlAk = NkN.k for A; ^ , 

AlAo = No{No - 1) < . (2.3) 



Then by ((221) and (Q one gets 

QIQ. = $^A*(A;)A(A:')A*A„ + 25^|A(A:)pA*A + |A(0)pASAo 

k,k'eA*, keA*, kytO 

= ^ E iyik)\{k')AlA,, + \%k')X{k)Al,A,) + 2Y,\m\'AlA,+ \\m 
^ I Yl ((^'^ + |A(A;)|)iVfc' + (iV_fc, + \X{k')\)N_k) + 2EiVfciV_fc + 

k,k'eA*, keA*, kjtO 

k^k', k^-k' 

= NkNk' + Y, NkN.k + N^+ J2\X{k)\Nk,. (2.4) 

k=^k' 

Using the inequahty 



fe.fe'eA*, fceA*, fc^O k,k'€A*, 

k^k' ki^k' , k^-k' 



NkN.k < \ {Nl + N\) , (2.5) 
we get 

NkN.,< (2.6) 
fee A*, fc^o fceA*, k^to 

Thus (EH) follows by (HH) and (HH)). □ 
We now use the inequality ()2.1|) in Lemma f2. II to prove superstability of the model (jl.8|l . 



Theorem 2.1 The Hamiltonian hl.^) is superstable: 

1 



- /^iVA >T^ + ^aiV' - (/i + i?)iVA (2.7) 
where R := Mu/2 and M is defined by (|i.^| ). 

Proof: From Lemma f2. II 

H^-liN^ > T^ + ^[v-u)Nl-{ii + R)N^ 

= T^ + ^aNl-{fx + R)N^. (2.8) 

Since we are assuming that a > 0, the estimate ()2.8p implies superstability, see [22] • D 
In the next two sections we develop the proofs for the variational formula for the pressure. 



3 The First Approximation 



Recall that the auxiliary Hamiltonians H/y{u,ri) and H^\q,v,rj) are source dependent with 
z/, ?7 G C, see p.22j) and p.2H|l . Since later we shall let v and 77 tend to zero, we can assume 
that |z/| < 1 and |?7| < 1. Because we are making the assumption on PBH p.8|l that m > 0, 
it follows from (IT^ that Hl{q) < 0. 

Let z/ G C and (j) := aigiy* \{k)) . Then from 

«±e-V,)(a,±e^Vj>0 

we get 

-\u\{Nk + iV-fc + |A(A;)|) < u\*{k)Al + u*\{k)Ak < \u\{Nu + N^u + |A(fc)|) . (3.1) 

Also 

Therefore, by Theorem 12. II one gets for \v\ < 1 and |?7| < 1, the estimate: 

H^{u,7]) ~ fiN^ > H^-J2{Nk + N^k + \m\)-N^-V-fiN^ 
> H^-{i2 + 3)N^-m^-V 

- T^ + ^(^Nl-{fi + 3 + R)N^-{M + 1)V . (3.2) 
Since Hl{q) < 0, we also have 

- T^ + ^(^Nl-{fi + 3 + R)N^-{M + l)V . (3.3) 

Proof of Theorem 11.21 : 

For simplicity we shall prove this theorem for u = 0. The proof for a general u follows 
through verbatim by translation for u ^ 0. Clearly since < 0, it follows from ()3.3|) that 
for any q we have for the pressure of the PBH (jl.22j) the estimate from below: 

Pa{v) > PA\q,j^ = o,?7) =PA\q,v)- 

Also for any q one obviously has the estimate from above: 

Pa(v) = PA\q,v) + (pa{i^,v) - PA\q,i^,v)) 



- iPA{i^,v) -Pa{v)) + {pA\q,i^,v)) -PA\q,v)) 
< suppi^^(g', 77) + (^Pa(z/, 77) - Pa ^(g, ly, 77) j 

- (Pa(i^,^) -Pa(^)) + sup (p^a\u, q',r]) -pA\q',v) 



and, therefore, we get 

snpp'l\q,T]) < j9a(?7) < suprf^(g,?7) + inf (pa(z/,?7) i/,?7)) 

q g 9 V / 

- {Pa{i^,v) -Pa{v)) + sup (pa ^ (g, z/, r/) -PA^(g,?7)) • (3.4) 

We shall prove in Lemma f3. II that, if z/a — as A | M*^, then 

liminf (pa(i^a, v) - PAiv)) = , (3.5) 

A 

and 



Next, with a particular choice of Ujy that tends to zero as A f R*^, we shall show also that 



limsup{sup(pi^(g, i^A,??) -pi {q,v))} = . 

A q 



limsup{inf(pA(i^A,?7) - Pa ^(g, i^a, ??))} = . 



(3.7) 



This last result (which is proved in Lemma fH.2j) is much harder and requires the arguments 
developed in [21]. Putting these together we get 



limpA(^7) = lim sup Pa (g, 77) 



that proves Theorem 11.21 . 

We now prove the two lemmas quoted earlier. 



(3.8) 

□ 



Lemma 3.1 



and 



liminf(j9A(«^A,?7) - PAiv)) = 

A 



limsup(pi^^(g, 1/a,??) -PA'iq^v)) = 

A 



Proof: Writing u = x + iy, using the convexity of the pressure and (jH.ll) we get 



d 

Pa{j^,v) - Pa{v) > ^(^^a(i^,?7) 



u=0 



+ y ( ^^a('^,^) 



u=0 



= y{^Ql + ^*QA)H,i,) 



fceA* 



> -^IH(2(iVAU( 



ttia ) > —K\u\ 



by (0|1 and Lemma EU Therefore if z/a ^ as A t M'', we get ^T^ : 

lim inf(pA(i^A,'7) - PAiv)) = . 

A 



(3.9) 
(3.10) 



(3.11) 



(3.12) 



Similarly one gets 

sup (p^A\q,'y,v) -PA\q,v)) < ^k|sup (^2 (A^A)//0(g,^,^) +"1^) < K\u\ , (3.13) 



by (ITU), (D and Lemma IbH Thus 

limsup{sup(pi^^(g, i^A,?7) -pA\q,v))} = , (3.14) 

A g 

that implies (jTTT]|l . □ 

Lemma 3.2 There exists a sequence {z^aIa ^^'^^ tends to as A f , such that 

limsupinf (pa(«^a,?7) -PA^(g, z^A,??) ) =0. (3.15) 

Proof: Using the Bogoliubov convexity inequality 

TT(A-B)e^ , ^ A . R Tr(A-fi)e^ 

^ ^ ^' — < In Tre^ - In Ttc^ < ^ ^ ' 3. 16 

Xre-B - - Tre^ ^ ^ 

and p.24|) we get the estimate 

< p,iu,r^) -p^^\q, u,r^) < {{Ql - Vq*)iQ, - ^g))^AM) ' 

Let 6Q^{u,r]) := - (<5a)//aM) 

AA(z/,r/) := (5g:(z/,r/) ^QaI^^, ^))haM) > 0- 

Then 

inf (pa{i^,v) -PA\q,'y,v)) < ^^a('^,^) • (3.17) 
We want to obtain an estimate for Aa(i^, ?]) in terms of u and V. 
Let 

where (■ , ■)h denotes the Bogoliubov-Duhamel inner product with respect to the Hamilto- 
nian H, see for example [21] or [23]. Using the Ginibre inequality (e.g. (2.10) in (23]) we 

get 

AA(z/,r/) <^{SQl{u,r]) 6Q^{iy,r]) + 6Ql{u,r]) 5Qa(^, ^))/fA{-,,) 
< D,iu,v) + i {PD,iu,v)}'/'{{[Q:, [H,iu,v) - /iiVA,QA]])HA(.,.)}'^' • 
We shall show in Appendix A that there is a real number C such that 

([Ql, [ffA(^,r/) -/^iVA,QA]])i^AM) ^ C'l^^/' . 

Thus 

A,iu,rj)<D,iu,rj) + {Cpy/'{V'/'D,{u,r^)Y^\ (3.19) 



From the definition of tfie Bogoliubov-Duliamel inner product we liave 



Here we consider tlie pressure Pa{i^,v) ^ function of two real variables, x = Dlcu and 
y = 3mi'. Since u > 0, then following the Approximating Hamiltonian Method for attractive 
interactions j21] we consider the integral 



d^dy—-pA{u, 7]) . 

2 OT^ 



With Z/+ := (5 + iy and z/_ := —S + iy, this integral is equal to 
Ia{S) = 

J [—5,6] 

Then by one gets 



-<5,<5] 



dx 



dx 



V 



lAm<- dy\(N, 



^5,6] 



where '■= X]fceA*(^fc + ^-k + Since by ()3.2|) and Lemma IB. 11 the expectation 

{Na/V) jj^f^^^-^ is bounded uniformly in z/ and in V , we obtain the estimate 



f d'^ 

I dxdy—pA{i^,v) 

J [-5,512 OX 



< 



2 
V 



-5,5] 



Similarly one gets the estimate 



-5,5]' 



92 

dx dy—pA{iy,r]) 



< 2C6. 



These give 



-5,5]2 



dxdyDA{y,r]) < CVb. 



(3.20) 



Since the integrand is continuous, by the integral mean-value theorem there exists a sequence 
{^a}a with \vf\ < (5 such that 



dxdyDA{v,r]) = (25)^L)a(i^a, 77) . 



-5,5] 2 

The last equation and inequality ()3.2Up imply that 

CV 



Da{va,ii) < 



A8 



which together with give the estimate 

Aa(i^a, r])<^ + 2^3/4^1/2 • 
Choosing 5 = (5a such that (5a 0, but V 6/^ ^ oo, we get 

lim:^ Aa(z/a,?7) = 0. 

By ()3.17p this completes the proof of the lemma. □ 
This proves the first approximation. In the next section we deal with the second one. 

4 The Second Approximation 

Note that from definitions ()1.2Hj) and p.27j) of the first and the second approximating Hamil- 
tonians, Hl^\q,i>,r]) and H^^\q, p,r]), respectively, it follows that 

H^^^ (g, iy = 0,r])- //f (g, P.v) = ^ v{N^ -pf>0. (4.1) 



Later in this section we shall show (see Lemma 14.11 and Remark 14.111 that 

(g,p,r/) <pf (|g|e'(-+2^),p,r/) =pf (|g|,p,r/). (4.2) 

In Lemma lOl we prove that for each g > there is a unique density p = Pa^Q^v) > 0) such 
that 

PA\q,pA{q,v),v) = inf Pa Hg, p, ??) . (4.3) 

p 

We can also show f Lemma 14. 3|) that there is at least one g = gA(?7) > 0, such that 

(2) ('2) (2) 

Pa {qA,PA{qA),v) = suppr(g,pA(g),?7) = supinfpV(g,P,?7) • (4.4) 

q q P 

For the sake of simplicity below we shall omit the variable r/, and we put 

Pa('?, v) ■= PA{q) and gA(?7) := gA . 
Finally, we shall show in Lemma [4.41 that if 77 7^ 0, then 

\im{p^^\q,,pM),v)-PA\qAe'^''-'''^'\v)} = 0. (4.5) 

A 

We start by proving Theorem 11.31 assuming the results of Lemmas 14 . 1 1 - which we prove 
later. 

Proof of Theorem 11.31 : 

We have to prove the limit p.28|) i.e. that 

(2) 

p{r]) :=limpA('7) = \imp\' {q^, PA{qA),r]) ■ (4-6) 

A A 



First, by ()4.H] and ()4.2|) we have for all values of the variational parameters g, p and the 
source parameter 77 that 

PA\q,v) := Pa H?, = 0, 77) <pA\q,p,r]) < PA\\q\, P,ti) . 

Therefore, 

PA\q,v) < infpr(|g|,p,?7) = p?\\q\,pA{\q\),v) 
p 

and thus by definition ()1.29|) we obtain 

supp'k\q,r]) < snpp^A\\q\,pA{\q\),v) = suppi^^(g, pA(g), 77) = Pa H^a, PA(gA), ??) • 

q q q>0 

This estimate implies that 

limsuppi^^(g,?7) < lipPA %a, PA(gA), ??) • (4.7) 

A q A 

On the other hand for all rj we obviously have 

snpp'-^\q,r])>p^l\qAe'^^+'^\v) = pf (^a, PaI^a), r/) (4.8) 

'pf (g-A, Pdq.),v) - P^A\qAe'^''^'^\v) 



Now the limit ()4.5p and the estimate ()4.8|) imply that 



limsuppi^^(g,?7) > lipPA ^(?a, PA(gA), ??). (4.9) 

A g A 



Taking into account ()4.7|) and ()4.9p we get 

limsuppi^^(g,?7) = limpA ^(^a, PA(gA), ??) . 

Combining this result with Theorem 11.21 we get ()4.6|) . i.e. the proof of Theorem 11.31 □ 
Now we return to proof of the lemmas quoted earlier. 

Lemma 4.1 Let the functions f and h and the spectral function E{k,q,p) be as defined in 
M.lb]) and il.l^) respectively. 

(i) If f{0,p) > u\q\ > 0, the pressure pA\q, p^rj) corresponding to H^\q, p,ri) is given by 

P^l\q,P,v) = --5^$^ln{l-exp(-/3E(M,p))}-i-^(i?(A;,g,p)-/(M) 

fcsA* ^ keA* 

i2f /(0,p)-|ng|cos(g-2V;) ) 1 , ,2 ^ 1 2 
+ 1 P(0,p)-n^|gP rr"^" +2"^ ' ^'-'^^ 

where 9 = arg q and tp = arg 7] . 

(a) If f{0,p) <u\q\, then pf^ {q^ P^v) infinite. 



Proof: (i) By ()1.16j) and ()1.27|1 we can write H^^\q, p,ri) — fiNj^ in the form 

H^A\q,P,v) - f^N^ = '^{f{k,p)alaf,-^{h{k,q)ala*_k + h*{k,q)a_i^a,^)} 

keA* 

- Vv{i^a* + v*a,) + VW{q,p), 

where ^ ^ 

W{q,p) = - - V • 

Let qX*{k) = \qX* {k)\e'^^''\ Then with = dkc''^^'''^/^ , for A; G A*, one gets 

HA\q,P,v) - = ^{f{k,p)dldk-^\h{k,q)\{ala*_i, + a^kak)} 

- Vv{r]e-''/% + r]*e''/^do)+VW{q,p), (4.11) 

where 9 = argg = 0(0). 

Note that if /(0,p) > u\q\ > 0, then fik,p) > \h{k,q)\ > for all k G A*, so that E{k,q,p) 
is well-defined and positive, see ()1.15|) . Let 

/(M , A ,2_i[ /(M n .412) 

Then the canonical Bogoliubov transformation: dj^ = x^a^. — yk<y^k^ gives 

(g, p, v) - P^A = X] E{k,q,p)alai, - W {^a^ + ^"0) 
fceA* 

+ ^Y1 iE{k,q,p) - fik,p)) + VWiq,p) , (4.13) 

fceA* 

where and , /c G A*, are boson creation and annihilation operators and 
We note that 

\e = \v 



2 _ I 12 /(O'P) - l^^l cos(e - 2^/^) 



E{0,q,p) 

From the diagonal form of Hj^\q, p,ri) — pNj^ in ()4.13p we get the pressure ()4.10|) . 
(ii) Now let /(0,p) < u\q\. Then the quadratic Hamiltonian ()4.1ip is not bounded from 
below. This means that the trace in p.l2|) is divergent and therefore the pressure pf^{q-, p, y]) 
is infinite. If f{0.,p) = u\q\, then by definitions p.l6|) and the conditions on X{k) at least 
the zero-mode term of the Hamiltonian ()4.1H1 is not positive. This again implies that the 
trace in expression ()1.12|1 diverges. □ 

Remark 4.1 From the explicit formula ^-lUj ) it follows that 

(g,P,^) <pf (|g|e^(-+2'^'),p,r/) =pf (|g|,p,r/). 
Recall that by M.2fJ\) and \4.10^ one gets for g > 

pf(g,P,^) = -^$^ln{l-exp(-/?E(M,p))}-^E(^(^^'^"^)-/(M) 

^ fceA* fceA* 



2 



bp 1 „ 1 . 



Lemma 4.2 Let r/ 7^ 0. Then there are numbers < pi{q,r]) < P2i(l,v) < ^> such that 
the infimum of p\ {q, p,ri) over p is attained in the interval {pi{q,v)y P2{q,v)) '^^^ ^/Pa(q') is 

/OX 

a value of p at which the infimum is attained, then dp\ {q,pp,{q),r]) / dp = 0. Moreover, if 
< go < 00, then 

inf {vpi{q,rf) - (// + uq)+) > and sup P2iq,v) < ^ > 
where s± := max(0, ±s) for s G M. 



Proof: By (j4.14p we have 



dpr , N _^ V / 1 /(M 1 / /(M \] 

+ vp. (4.15) 



f |?7| 



From (prTH|l we sret 



(/(0,p) - ng)^ 

-(g, P, V) < - frfn \ ^ + ^P 



Let X := vp — {p, + uq)+. Using the identity fi + uq = (/i + Mg)+ — (/i + wg)- we obtain 

As X — * 0, the right-hand side of the last inequahty becomes negative. Therefore, there exists 
6{q,rj) > such that the infimum of p^l\q, p, r]) over p cannot be achieved if vp — (fi + uq)^ < 
d{q,ri), i.e. p < pi(g,r/) := ((p + Mg)+ + 6{q,ri))/v. 

It is clear that if < go < cxd, then infq<gQ S{q,v) > 0- 

Suppose now that p > pi{q,r]) and take wp > max(2p, 2g + 2). Then for G A* one has 
E{k,q,p) > max{e{k) , 1) . Therefore, using 

- Eik,q,p) - E{k,q,p) " ^' ^ 

we obtain the estimate 

.(2) 



^^W,P:^ = -77 i TT^TTTl 7 + o coth A:,g,p — --1 

9p [ exp{(3E{k,q,p)) - 1 2 2 \E{k,q,p) J 

vp (4.16) 



fceA* 
t>|?7| 



(/(0,p) - ng)2 



V exp[/?max(e(A:),l)]-l 2V ^ ^' 5(g,r/)^ 



wp 



Making use of ()1.4j) . this implies that there exists a volume Vq independent of g and p, and 
K{q,r]) > such that if > Vq, then 



-^iq,p,v) > -Kiq,v) +vp, 



and therefore, if p is large enough, then — — (g, p, 77) > 0. As a consequence, there is P2i(l,v) 

op 

such that the infimum of Pa {(l,P,v) is attained in the interval {pi{q,ri), p2{q,ri)). If Pa(q') is 

a value of p at which the infimum is attained, then — — (q, PAi(l),v) = 0- 

Let < go < 00. Then one can see that supg<qjj K{q,rj) < 00, and therefore we get 
supQ<qoP2(g,^/) < 00. □ 

Lemma 4.3 Let r] 7^ 0. Then there is qo{ri) < 00 such that the supremum of p\ [q^ p ^^{q) , r]) 
with respect to q is attained in the interval (0,go('7)) for all A and if q^ is a maximizer of 
PA\q,pA{q),v), then 



dq 

There exists co{r]) such that for all A 



(gA,PA(gA),^) = 0. 



/(0,PA(gA)) - uqA > 00(77) , 
if (Ja is a maximizer of p^^\q, pj^{q),r]) . 

(2) 

Proof: Recall that v — u := a > 0. Differentiating pX (g, p, 77) we get 

^(g,p,7/) = -^g|A(A;)p| g^p(^^(^^^^^)) _ ^ ^(k^^) + 2E{k,q,p)] 



|2 



U\T}\ 



By Lemma f4. 21 we have 



dpf^, _ , . . _ _ .dp^{q) dpf^ _ 

-^{q,pA{q),V) = -g^il^pAiqj^V) + -g^iQ.pAiqj^V) = -Q^il^pAiqj^V) , 

since dp^^\q, p/^{q),ri)/dp = 0. Therefore, we can also write 

dp^^ dp^^ dp^^ 

-^(9, pA{q),v) = ^^(^' PA(g), + ^^(?' pA{q),v) ■ (4-18) 
Insertion of ()4.15j) and ()4.17j) into the identity ()4.18j) gives 

dp^A\ , 1 r 1 vf{kM(i))-u'<i\m? 



1 

iq,pA{q),v) = -y^{ 



dq ' ' V exp{(3E{k,q,p^{q))} - 1 E{k,q,p^{q)) 

1 f vf{kMq))-u\\x{k)\' 



{fioMq))-^qy 



+ vp^{q)-uq. (4.19) 



Then, since f{k,p) > uq\X{k)\ > uq\X{k)\'^, f{k,p) > E{k,q,pj^{q)) and a > 0, by ^UH^ we 
get the estimate 

dpA^ , ^ ^ 1 u'^q\X(k)\'^ akP , . , . 

E{k,q,pA{q)) (/(0,pA(g)) - uq^ 

Now we have 

E\k,q,p) = {f{k,p)-uq\X{k)\){f{k,p) + uq\X{k)\) 
= (e{k) + {fiO,p)-uq} + uq{l-\X{k)\}) 

x{e{k) + {f{0,p)-uq} + uq{l + \X{k)\}) 
> {f{0,p)-uq)uq . (4.21) 

Therefore, by dOj), dm and Kld^ . KTOi we obtain 

—^{q,pA{q),V) < oT//f^n - f \\ \m~ fffn - f \\ T^+f{^,pA{q))-uq+p. 4.22 

dq 2V{f{0,pA{q)) -uq)y^ {f{0,pA{q))-uqy 

Let aA^q) '■= (/(0,Pa(q')) — uq) (max(l, g))^/^. Then the inequahty ()4.22j) gives 

dp?\ , ^ (mm(l,g))'/3 [ cMm'/' a|^P 3/2, , . „„, 

dp?^ 

Therefore, there exists Cofr/) such that if g > 1 and crJq) < Coir]), then — ; — {q^ pAq),^) < 

dq 

for all A. Thus for all A the supremum of Pa^^ {q, PA{q)-,v) o^^r q cannot be attained in the 
domain defined by the condition 0"A(g) < Co(?7). 

Now assume that q > 1 and Cj^^q) > Co(?7). Then, using again ()4.2H) . we obtain from ()4.17j) 
the estimate 

dp^A^ r 1 q^^ 

-iq,pAiq),v) < K 



dq ^--^^^^"^ - l(/(o,pA(g))-«g) {f{^Mq))-uqYi^ 

2 



Since the right-hand side of ()4.24j) becomes negative for large g, there is gol??) < oo such 

('2'] 

that the supremum of Pa {q, PA{q),v) with respect to g is attained in g < go(?7) for all A. 
Note that from ()4.17j) we see that if gA is a maximizer oi px {q, PA{q),v)^ then gA 7^ 0, and 
therefore combining this with the last statement we can deduce that 

^(gA,PA(gA),^) = o. (4.25) 

Putting 00(77) = Co(?7)/{max(l, go(?7))}^/^ finishes the proof. □ 



Lemma 4.4 // 77 7^ 0, then 

lim{pi')(g^,p4gj,r?) - rf^lgAe'^^+'^^U)} = . (4.26) 



A 



Proof: By Bogoliubov's inequality ()3.16|) one gets 
Let 5Nj^ := — ^Pa(^a) and 

■■= (^^^A><)(,.eu....),,,(,,)„) • (4.28) 

Then ^717^ implies 

0<p'f\q.,pdqA),v) -p'A\q.e'^''^'^'\v) < ^^a(^) • 
We want to obtain an estimate for Aa(?7) in terms of V. To this end we introduce 

DM = (4.29) 
and calculate the derivatives 

.(2) 



dp':' , . ^ 1 [ 1 /(M 1 / /(fc,p) 

P, ^ \ exp{l3E{k,q,p)) - 1 E(A;,g,p) 2 U(^,?,P) 



(4.30) 



(/(0,p)-«g)2' 
(?,P,^) = -^'H7^(9'P'^)-P ' (4.31) 



dp \ 9/i 



^2 (2) 

-(9,p,^) 



T/ 



/3exp(/5E(M,p)) /2(A;,p) , 1 exp(/5E(M,p)) + 1 «V|A(fc)| 



^ fct^' (exp(/3E(A;,g,p)) - 1)' E\k,q,p) 2 exp(/5E(/c,g,p)) - 1 E\k,q,p) 



(4.32) 



(/(0,p)-ng)3 

From (lO^ . using e^/(e^-l) < 2(1 + 1/x) for a; > and P{k,p) = E{k,q,pf + u'^q^\\{k)\'^, 
we get the estimate 

^^^'^'''^ " ^£ (e^p(/5^(fc>?>p)) - 1) v^W^) 



4e[ 



2PE{k,q,p) + 3 1 



y I {exp{PE{k,q,p)) - 1) 2E\k,q,p) 

2|r/| 



(/(0,p) - ug) 



u'q'\m? 

(4.33) 



The second sum in fl4.HH|l is bounded from above by 



V i^^\E^(k,q,p) E'(k,q,p)J \(f(0,p)-ugf {f{0,p)-uq) 

and the first sum (using E^{k,q,p) > e{k){e{k) — /i)) by 



01 



V ] f^, \Eik,q,p) E\k,q,p)J (exp(/56(A;)/2) - 1) 

e{fe)<l + 4|M|/3 £(fe)>l+4|M|/3 

- ((/(0,p)-ng) + (/(0,p) - + ^ 
Consequently 

where Co(?7) and q'o('7) are as in Lemma f4. 31 

By LemmaOwe have ^(g., p.(g.), = 0. Then from ^ one gets that 

PaI^a) = %^(gA,PA(gA),?/) = ^^(gAc'^^+^'^^PAlgA),^) = 



(4.34) 



<-) VIA) fii VlAy ) ■/ y r% VIA" ^fAVIAyj'/y \ / ' 

\ ^ /<\5a,Pa(5a),.) 

and therefore by ()4.29p 

It then follows from ()4.34j) that 

lim^ = 0. (4.35) 
Now Ginibre's inequality for ()4.28p and ()4.29p . cf. Section 3, gives 

^Aiv) < DM + (4.36) 
-P'/'{dm] ([A^a, [i/P(gAe^(^+^'^),p,(g,),r/),iV,]]) . 

2 I \ / ^^r(9Ae'(-+2'/'),PA(gA),r;) J 

Note that here 

([iV„ [H^^\q,p,r^),N,]])^^^^^^^^^ = 2u {q*Q, + gQD^c^)^^^^^^) + V^(r/a* + r/*ao)^(.)(^_^_^) • 
By differentiating the pressure we find that 

2V ( dpi^ dpi^ 



so that if we define q := \q\e^'^'^^'^'^\ then we get 

Hi (I'P'V) u \ ' ' d\q 



4,V ( dp^f' 

u{q*Q^ + qQl) rrm(,,„. = 2u\q\'V + — \q\——{\q\,p,r]) 



An exphcit calculation gives 



drj drj' 

f{0,p) - u\q\ cos{e - 2tp) 



/2(0,p)-n2|g|2 



and so 





1^1 


2 


(/(0,P) 


— u\ 





('K + 'r«o)„«.«,„„ = [jf^jj^^ \ . (4.37) 

Therefore, if dp^^\\q\, p.rj) /d\q\ = 0, then 

([iV„ [i/f (g,p,r/),iV,]]) =2V (2u\q\' + 

Thus 

From ()4.35|1 . ()4.36|1 and the last estimate we then see that 

completing the proof. □ 

Now we prove that the order of the thermodynamic limit and taking the infimum and 
supremum in ()4.2fij) can be reversed. 

Proof of Theorem 11.41: 



We know from Lemma 14.31 that there is q'o(^) < oOj independent of A, such that for 
large A, the maximizer ^ [0, q'o('7)]- Then it follows from Lemma 14.21 that ^^{ri) : = 
infqG[o,go(»;)] vpi{q,ri) - (/i + uq)+ > and po2(r/) := supgg[o,go(»;)] P2(g,r/) < oo. Thus pA{q) is 
in [0,po2{ri)] and vpA^q) — (/i + > 5o(^)- Let C be the compact set 

Gri ■■= {(g,p) I < g < qo{r]), [{p + uq)+ + 6o{ri)]/v <p< po2(^)} • 

Then (gA,PA(^A)) G C,,. Therefore, there is a sequence A„ such that (^a„, Pa„(^a„) converges 
to some point (g, p) in G^. 

The derivatives of Pa^^ {q-, p,ri) are uniformly bounded on and therefore as A f W^, 
p^A\q,p,v) converges uniformly to p^'^\q,p,rj) on G^. Thus 

limpA(^) = lim pS(^a„, PaJ^aJ, ^) =P^'^\q,P,v) ■ (4-38) 

A n— >oo 



By repeating the arguments of Lemmas 14.21 and 14. HI and by replacing (for V oo) the sums 
over k by integrals, we see that if q maximizer of infp.o-(q,p)>op''^H?, P, "'?) with respect to q, 
then < q < qo{v) if p{q) is a minimizer of p^'^\q, p,r]), then {q,p{q)) is in G^. Thus 

sup inf P^'^\(li p,v) = sup inf p^'^\q, p,T]) . 
q>o p-<^{q,p)>o gep.goW] {p^I'J.p)^^,,} 

Since 

(^An . PAn (^A„ ) , ^) < (^A„ , P, ^) 

for p such that {qA„,p) £ G^, we get also that 

for p such that {q,p) E Gr,- That is 

P^^H^, P, ^) = , inf P^^\q,P,v) ■ (4-39) 

Similarly, for all g > we have 

Pa! {qA„ , Pa„ (gA„ ) , ^) > (g, Pa„ (g) , ^) • (4.40) 

If < g < ^0(^7)! then (g, PA„(g)) G and therefore Pa„{<i) has a convergent subsequence 
Pa„,, (g) converging to some p, where (g, p) G G^. Taking the limit in ()4.4Up we obtain 



P 



(g, p, ^) > p^^^ (g, p,v)> , , inf ^ , p^^^ (g, p, r/) • 



{p:(g,p)GG^} 

Thus, by 

P^^^ (g, P, ^) = , , inf ^ ^ P^^^ (g, P, ^) > ^ ^ inf p^^^ (g, p, r/) 

for all q G [0, 90(^7)] • Therefore 

p^'^\q, p,r]) = sup inf ^'•^^(g, p, 77) = sup inf p^'^\q, p^r]). (4-41) 

gep.goW] {P-(<l'P)^G^} q>o p:cr(q,p)>0 

Combining the relation 1)4.411) with ()4.38p we prove the theorem and obtain an explicit 
formula for the limiting value of the pressure. □ 



Remark 4.2 By the definition of p^iq), q^ (Lemma \4-S\\4-3}) and by l{4-38 ) we also get that 
for rj ^ 

lim p2 (gA„ , Pa„ (gA„ ) , ^) = lim sup inf {q, p, r]) = sup inf p^^^ (g, p, r]) 

n— >oo n^oo g>Q p>0 |j>o p>0 

= p^'Hq,p,v), (4.42) 
where (of. I{4.14\ )) 

p^'\q,p,v) = - [ yf^(^ln[l-exp(-/3E(M,p))] + ^(i5(M,p)-/(M) 

+ TTTTT W + W^ (4-43) 

/(0,p) - ng 2 2 

and g, p satisfy the equations 

^{q,p,v) = , ^(g,p,r/) = 0. (4.44) 



We now show that the zero-mode ?7-source term can be switched off. 



Lemma 4.5 Thermodynamic limit of the pressure is equal to 

p := hm »A = hm hm pAn) . 



Proof: By Bogohubov's convexity inequahty ()3.16p one gets 



that imphes 

< \p.-pM\ < %\ I < ^ («S«o)i(,) < 4= • (4-45) 



From Lemma [B. II and ()3.2|) we see that for |?7| < 1, 



Ha{v) 



V 



where Ki is independent of r]. Thus the right-hand side of ()4.45j) tends to zero as t] tends to 
zero. □ 

Finally we prove that the order of the limit ?7 — and taking the infimum and supremum 
in ()4.4H1 can be reversed. 

Lemma 4.6 

limsup inf p*^^-* (g, p, 77) = sup inf p^'^\q,p) , (4.46) 

q>0 P '■ o"('?.P)>0 g>0 P : o-{q,p)>0 

where p^'^\q^ p) := p^'^\q, p,0) is defined in \1.14}) . 

Proof: Let Priiq) be such that 

inf p^^\q,p,v) = p''^\<l^Pvi(l)^v) , 

and be such that 

supp(^^(g,p^(g),r/) = p^^\qr„ pr,{qv) , v) ■ 

q>0 

Let 

Go :={(g,p) \ q>0,a{q,p)>0} . 

By arguments similar to the above (see proof of Theorem II. 4p we can show that these exist 
and that {q-q, Pr^iqr^)) £ Go- We shall need the following derivative of ()4.43|) : 

dp^'\ . f d'^k f 1 /(M , 1 f f{k,p) 

v\ri\ 



dp " U C^T^r I exp{(3E{k,q,p))-l E{k,q,p) 2 \E{k,q,p) 

2 

vp . (4.47) 



(/(0,p) - uqy 



Moreover, in the same way as in ()4.17j] . ()4.19|1 we also obtain: 



dq 



.(27r)-''^' { exp{f3E{k,q,p,,{q)))-l E{k,q,pr,{q)) 2E{k,q,pr,{q)) 



(/(0,P,(g))-wg)^ 

and for any number t 



uq, (4.48) 



dp^'\ , f d-^k j 1 tvfik,p,iq))-u'q\Xik)\'^ 



dq ' ' JA'^ttY 1 exp{(3E{k,q,Prj{q))) - 1 E{k,q,pr,{q)) 

1 ftvfik,p,iq))-u\\Xik)\^ 



2 V ^(A;,g,p,(g)) 
a|r/|2 



tv 



{f{^Mq))-uqY 

As in fl4.24|) . from (j4.48|l we get the estimate 



+ tvp^{q) -uq. (4.49) 



^(g,P.(g),^) < ^ I (/(0,p,(g)) - ng) + (/(0,p,(g)) - ng)V2 

^ (/(0.p,M)-„.)^ -"^- 

Therefore, if /(0,p^(gr/)) — wg^ > 1, then by the definition of and by ()4.5()j) we obtain 

Since the right-hand side of the last inequality must be non- negative, then 

/(0,p-,(g,))-ng,<^iii^t^. 
Similarly, if f{0,pr^{qr^)) — uqn < 1, then 

The right-hand side of the last inequality must be positive and thus 



f{0,pr,{qrj))~uq^< 

U^i'^qri 



Therefore, either 



1 < /(0,p,(g,))-ng, < ^ or < /(0,p,(g,))-ng, < min 1, ' \^ % ' 



Thus the only way that PrjiQrf)) can escape to infinity as r/ — ^ is, if either Pr]{%) oo 
and q^i 0, or if Prjigr^) ^ oo, ^ cxo and /(0,p^(g^)) — uq^ — > 0. Now, if p — > cxd and 
q —>■ 0, the right-hand side of ()4.47|) tends to +oo. Therefore the case pr^{qr^) — oo and 
qn — > 0, is not possible. 

Suppose now that Pr]{qri) — > cxd, — > oo and /(0,p,,(g,,)) — uq^j — > 0. From ()4.49|1 with 

t = u/v we get 



dp^'^\_ , , ||A||n _ _ ||A||m ^^/./^-/-xx - -X 

= ~dq~^^'" P'yl^r,), ^) < + up^[qn)) - uqr, = + - (/(0,p^(g^)) - uqr, + p - aqr,) . 

This contradicts our supposition and therefore p-qiq-q) and must remain finite. 
As in fRT^ and (jOSl), from (jOn|l with t = 1, we get 

^ dp(^), - /- N N 1 / 11-^11 ^^^^^y^ 

= -3— (g,y,Pr,(g,;),^) < 



dq " (/(0,p,fe)) - uq.yl^ \ 2 (/(0,p,fe)) - uq.fl^ 

/(0,Pr,(gr,)) - Uq^ + P. 



Therefore, since the right-hand side must be positive, the term 



{f{0,pM))-uq,r/' 
must remain bounded when /(0,p^(gr;)) — uqq 0. 

Summarizing we see that {%iPr]{%)) must remain in a hounded subset of Gq and 
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lim , _ ^' — = . (4.51) 

''^o {f{0,Pri{qn) - uq-n) 

Since {qri, Priiqrt)) remains in a bounded subset of Go, there exists a sequence r/n — ^ such 
that {qrin, PriniQr)^)) converges to (q, p) G Gq; where Gq is the closure of Gq. Now p^^\q, p) is 
continuous on Gq. Thus by ()4.51|) we obtain 

limpA = limp(2),^^ pj^j^^^) 

A n^oo 

1^12 

= lim (g^^ ^ (g^ J ) + lim 



= p^^\q,p) ■ 

Now for p such that (g, p) G Gq, for large n we have {q^^, p) G Gq. Therefore, for large n we 
get 

and letting n ^ 00, we obtain for p such that (g, p) G Go, the estimate 

P^'^\q,p) < P^'^\q,p) ■ 

That is 

p^^\q,p) = , , inf P^^^(g,p). 



Similarly, for all g > we have 

From ()4.47j) . we see that for each q > 0, both Pri{q) and \'r]\'^/ {f{0,Pri{q) — uq)'^ remain bounded 
as ?7 — > 0. Let {Prin^il)}^ >^ be a convergent subsequence of {Pr]n{(l)}n>i converging to p 
say, where (g, p) G Gq. By letting r ^ oo we then have 

P^^\q,p) >P^'^\q,p) > ^ inf P^^\q, p) ■ 



Therefore 



p^'^\q,p)= inf p^'^\q,p)> inf P^'^\q,p), 
{p:{q,p)eGo} {p:(g,p)eGo} 



for all g > 0, and thus we get the relation 

p*-^-* (g, p) = sup inf p^^'' (g, p) = sup inf p^'^\q,p) 

q>0 {P-iQ>P)&Go} q>0 P-o{q,p)>0 

proving the theorem. □ 

Combining Theorem II. 41 Lemma and Lemma we get the first part of our main result. 
Theorem O (11201) • 

The second part we shall consider in the next section. 



5 Discussion 



Let us put in Hamiltonian (jl.22j) the source equal i/ = and suppose that i] 0. Then 
the corresponding Euler- Lagrange equations, obtained by the condition that the derivatives 
()4.47|) and ()4.48|) are equal to zero, take the form 

<!=-[ 7^#^ cothV(M,p) + 77^^^- (5.2) 
2 U (2vr)- E{k,q,p) 2^ ^ '^'^^ (/(0,p) - uq)^ ^ ^ 

We shall now discuss some of the consequences of these equation in relation to the existence 
of Bose-Einstein condensation (BEC) in the model ()1.22p . 

(a) The solution {priiP, p), q-qiP, p)) of the equations ()5.2|1 always exist and is a smooth 

function of f3, p and 77, for rj ^ 0. Moreover, we can identify it with the Gibbs expectations 
of the corresponding observables. Since the pressure Pa^i' = 0,r]) is a convex function of p 
and of u, then by the Griffiths lemma, see e.g. jlj, the corresponding derivatives converges 
in the thermodynamic limit to derivatives of the limiting pressure fjl.30|) . Differentiating 
fll.H()|l with respect to p and u and comparing these derivatives with the solutions of (j^.lj) 
and ()5.2j) . we get 



(b) Similarly we can show that the zero-mode BEC for rj ^ is given by 

Po(r7):=lim/^\ = .„.^Jf _ ■ (5.3) 

To obtain this result let us make a global gauge transformation U^p = e**^^* of the Hamiltonian 
Hf,{n, z/ = 0, r/) = Hj^{iy = 0, 77) — /iiVA, see p.22|l . with cp = argr]. Then : 



Hj,{jj,, 0, 77) = UpHj,{jj,, 0, ?7)t/* = Hp,- /iiVA - Vl/|?7| (flo + ^0) • 

From 

= ([^A(/i,0,r/),iVA])j^^(^o,^) = VV\Tl\{al - 00)^^(^,0,,,) 

and 

< ([iV,J^,(/i,0,r/),A^J])^^(^^o,.) = v^I^I(«o + «o)h.(mA.) 

we obtain 

(«S)HA(M,Ofl) = («o)ha(^,o,„) > 0- (5-4) 
Let (5Ao := {al + Oq) - (fi* + ao)^^(^ o ,,). Then 



(^^S, ^^o)^,(,,o,.) > 0, (5.5) 



where (■ , OffACM-o,??) denotes the Bogoliubov-Duhamel inner product with respect to the 
Hamiltonian Hf^{ix,v = 0,7]). Hence, the convexity ()5.5p and convergence of the pressure 
Pa{v) (see Theorem 11.41 and Remark I4.2jl imply by the Griffiths lemma the convergence of 
the first derivatives to the derivative of the limiting pressure : 

r ^Pa(^) 2\ri\ 

hm = lim —= {ar. + an) fv^ f,^ n ni = ./^ - n — , (5.6) 



see fll.30|) . ()1.3ip and ()4.37|) . Therefore, by ()5.4p . ()5.6p . and returning back to original 
zero-mode operators, we obtain 



lim ( ) = _ ' lim ( — - ) = _ ' (5 7) 

So, by ()5.7|) we conclude that the r] - source in Hamiltonian p.22p breaks the zero-mode gauge 
invariance creating a zero-mode macroscopic occupation with the particle density estimated 
from below by the Cauchy-Schwarz inequality: 

lim/4^\ > lim/4\ /^\ (5.8) 

^ \ V / Ha{o,v) ^ \VV/ha{o,v) \vy/ hm 



\V 



2 



(/(0,Pr,) - Uqri) 



To prove that in fact there is an equality in (j5.8|l . we consider p^iv^ ^) ^^e pressure with e(0) 
replaced by e(0) — s with s positive and again use its convexity with respect to s. Then 
Griffiths lemma and that fact that f{0,pr^) —uq^j > 0, as soon as 77 7^ 0, imply, see ()4.14|) and 
dOHl): 



\ =\im ( ^E^^i^Zl^^ < f M^Zl^^ 



Here we have used the fact that the s-dependence of p(?7, s) is only through the last term in 
(HSH). 



(c) In the limit rj ^ equations (jS.lj) and ()5.2j) coincide with equations (3.7) and (3.8) or 
(3.10) and (3.11) in [13^. There the amount of the generalized condensate density is denoted 
there by mo(/3, /i). By inspection this coincides with the limit of po{v) in ()'^.3|1 as r/ ^ 0: 

mo{(3,^i) = linipo(^) • 

In we found that for mo to be non-zero, fi must be greater than a certain critical value 
of chemical potential Hc{f^,u,v) . For u = 0, this critical chemical potential coincides with 
the one for the Mean-Field boson gas ()1.10|) . namely ^c{P,u = 0,v) = vpdP), where pdP) 
is the critical density for the Perfect Bose-gas, see e.g. 0. 

(d) It was shown in [TH] that the phase diagram is quite complicated. Subject to these 
Euler-Lagrange equations the expressions for the pressure given in jTHj equation (2.11) and 
at the top of page 438, are the same as p^'^\q, p) in ()1.14|) . (We warn the reader that in these 
equations for the pressure in jTSj there is a misprint and a term is missing.) There we were 
able to solve the problem only for some values of u and f , see Fig. 2 in ^HI- For example 
()5.2|1 shows that for n > [attraction in the BCS part of the PBH ()1.8|) ) the existence of 
the generalized Bose condensate mo 7^ causes an abnormal boson pairing: 

limlim ^ {Ql + Qa)ha(o,^) = Jim ^r,(/5, /i) ^ . (5.9) 

This is because, for n > 0, equation ()5.2|1 cannot have the trivial solution g,, = when the 
generalized condensate 

mo(/3, /i) = lim — ^ ^ . (5.10) 

(/(0,p,,) - uq^f 

Note that on the other hand the equations 1)5.11) and ()5.2|) allow the possibility that mo = 
without lim^^o % = 0. This "two-stage" condensation is possible only when u > and it is 
similar to that discussed in [T!?j . 

(e) As in we interpret the spectrum ()1.15|) of the effective Hamiltonian 

£e^citik) := lim E{k,q p ) , (5.11) 

as the spectrum of excitations for the PBH ()1.8)1 . Our analysis of the Euler-Lagrange equa- 
tions ()5.1|) . ()5.2|) (as well as ()5.13p . ()5.14l below) shows that there no gap in this spectrum 
as soon as there is the Bose condensation ()5.10p : 

limeexcit(fc) = limlim(e(/c) -fi + p^- \uq^\*{k)\) =0. (5.12) 
This conclusion is again in agreement with [T!^ . 

(f) The case of repulsion {u < 0) in the BCS part of the PBH (|1.8p is quite different. In 
this case the pressure coincides with the mean-field one {u = 0) and we always have for the 
boson paring: lim^_>o qr^iP, p) = 0. The first property was derived in great generality in [T!^ . 



To make a contact with the variational principle proved in this paper, let us change notation 
and replace u by —w, with w > 0. The Euler-Lagrange equations, ()5.1|) and ()5.2|) . become 

I coth^/3E(fc,g,p)-l| + ,,,, y ^, (5.13) 





/ a k 








{-w)q 1 


Q = 


2 k 



r ^ '^'"^ j (/(0,p)+.^g) 
' ^ ^' coth + 7777-^^ rrj • (5.14) 



. (27r)- E(M,P) 2" ^ '"'^^ (/(O, p) + i/;g)2 " 

Since the solutions p,,(/5, /i), g^(/3,/i) of equations (j5.13j) . (j5.14j) must satisfy the condition 
cr(g^,p,,) > 0, one gets by ()1.19j) the estimate 

/(O, p^) + w% > 2wq^ . (5.15) 

Note that the first term in the right-hand side of ()5.14j) is negative. Therefore, by ()5.15|1 we 
obtain 



\v\ \v\ - \v 

g^(/3,p)<— — — — TT^ < TTT— ^ (lr,{P,fi)< 



2/3 



(/(0,p)+^g,(/5,p))2 {2wq,{P,fi)y ^^^""^^ (2.^)2/3- 

This implies that in the limit 77 — > the equation ()5.14j) may have only a trivial solution: 

limg,(/5,p) = , (5.16) 
77— >o 

and 
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lim— 4^ ^77 = 0. (5.17) 

Let Pc(/9) be the critical density for the Perfect Bose Gas: w = f = 0, see p.9p or ()1.10|) . 

dJ'k 1 



Pc(/9) : = 



For p < vpc{l3), limits ()5.1(ij) and ()5.17j) imply that as 77 the solution of equation ()5.13j) 
tends to p(/3, p) the solution of the corresponding equation for the Mean-Field model ()1.1U|) : 



P 



d^k 



. {2tiY e^(^('=)-^+''/') - 1 



and the pressure 



i5'^(/?,p):=lim inf ird p^^\q, p, 7^) = inf p(2)^o, p, 0) = p(2)(0, p(/5, p)) 

ri-*0 p :a{q,p)>0 q>0 p:fT(0,p)>0 

coincides with the mean-field pressure, see ()1.14|) and jH]. On the other hand, if p > ji/v, 
then from ()5.13j) we obtain for any e > and 77 is sufficiently small 

^ < P7?(/5,p) = Pc{^)+e , 

giving a contradiction for p > vpc{f3). This means that in this case equations ()5.1H|) and 
()5.14|1 are inconsistent and the minimum point must lie on the boundary of the allowed 
range on the p-q plane. This boundary consists of the two lines q = and p = {fi + wq)/v. 
Minimizing the pressure on the first line is equivalent to solving the variational problem in 
the mean-field case. This was done in j5] where one sees that the minimum is attained at 



a point which tends to p = fi/v as rj ^ 0. On the other boundary p = (fi + wq)/v similar 
calculations show that the minimizer also tends to [p = fi/v,q = 0). Thus the pressure again 
coincides with the with the mean-field pressure. 

This proves the second part of our main result for repulsive BCS interaction in the PBH, 
Theorem EH |T^ . 

We end with the following remark concerning BEG in the PBH model. Though the pressure 
of the model with the PB Hamiltonian for w > coincides with the one for w = 0, it is an 
open question wether these models coincide completely. As has been shown in [221- [SI] ^ 
similar type of diagonal quadratic repulsion is able to change the type of Bose condensation, 
from condensation in the zero mode (type I) to generalized van den Berg-Lewis-Pule con- 
densation jnni out of the zero mode without altering the pressure. Therefore, the analysis 
of the Bose condensate structure in the PBH model requires a more detailed study of the 
corresponding quantum Gibbs states. This is beyond the scope of the present paper. 
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Appendix A: Commutators 
By (Hini) and (fT^ we have 

[H,{u,v)-fiN,,Q,] = (-2) Y,m-f,)X{k)A, + y Y.\X{k)\'(N, + ^) Q, 

keA* keA* ^ 

- ^{N.Q, + QM+Au Y,\m\'(Nk + l)-2 

keA* ^ ^ 



Vrja, 



' 



and 



+ 



fceA* 

y I E \X{k)\'y{k)Al Q, + 2 IM^)!' (iV. + I) 
I fceA* UeA* ^ ^ 

QlQ^ + 2Y\m\'(Nk + l) (iV^ + l) 
fceA* ^ ^ 

8z/ E |A(A;)pA*(A;)A* +4v^r/A(0)a*. (A.l 



V 



fceA* 



Using p.5|) and ()1.6|) we see that the first term in ()A.1|) is bounded by 

8(ca + \p\){N^) + 4nA + 4|/i|mA , 
where ( ■ ) := ( ■ )hj^{u,ti)- Recall that Lemma ED gives 

QIQ^ < Nl + MVN^ 



and as in (jSISI) we get Af^Al < Nf^N_^ + 3{N^ + A^_^.) + 2. Using these we obtain 

fcGA* fcGA* 

VfceA* / VfeeA* 

. 1/2 

< (iV,2 ^ MV^iV,)!/^/^ ( \m\{N,N_, + 3{N, + N_,) + 2) ' 

VfceA* 

< {Nl + MVN,)^'^xnT {{Nl + 6iV, + 2m,))'/' , 
and independently we have 

TTIa 



fceA* ^ ^ fceA* ^ ^ 

and 



5^ |A(fc)P (iV. + ^) (iV. + 1) < E 1^(^)1 (^^ + ^) + 1) ^ (^^ 

fcGA* ^ ^ fceA* ^ ^ 

which gives estimates for the second and the third terms in ()A.1|) . We now bound the 
penultimate term in ()A.1|) . 

fceA* fceA* fceA* 

VfceA* / VfceA* / 

Finally for the last term we have 

Putting these bounds together we get 



y{Nl + MVN^Y'^xnT {{Nl + 6iV. + 2m.))'/' 
8m mA\2 Av 



+y[^- + ^) +y{Nl + MVN,) 

+8 {N,)'/' {{N,) + m,f' + 32VV {N,)'/' . 
for \u\ < 1 and |?7| < 1. From Lemma fB. II and ()3.2j) we see that for < 1 and |?7| < 1, 

(^) and (^\ <i^„ (A.2) 

where Ki and K2 are independent of 1^,77. Thus 

([Q:,[i/.(z/,r/)-/.iV„QJ])^^(,_^) <C\/='/2 (A.3) 

for some number C. 



Appendix B: Bounds 



Lemma B.l If a Hamiltonian satisfies the condition 

Ha>Ta + ^^Nl - 6Na - aV (B.l) 

with 7 > then there exist constants Ki and K2, depending only on 7, 6, a and /x but not 
on A, such that 

^) < K, (B.2) 



V 

\ 



and 

<K2. (B.3) 

Ha 



'A 



Proof: Let Pa(/^) be the pressure for Ha, then 

'iVA" 



. . / < Pa(/^ + 1) - PA(/i) < Pa(/^ + 1) < ^1, 

^ / Ha 

where Ki is independent of A by (jB.ljl . Also for AG [0,7) let 

Ha{X) :=Ha-^XNI 
and let Pa(/^; X) be the corresponding pressure. Then 

T7I ) < -{pA(/i,7/2) -pA(/i)} < -PA(/i,7/2) < 
V / Ha ^ 7 



where is independent of A, again by ()B.1|) . □ 

Note that by Theorem 12. II the Hamiltonians (jl.Sj) and p.22j) verify the condition (jB.lll . see 
estimate ()H.2|) . 
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